attach(coking)
anova(lm(time ~ width * temp, data = coking))
tapply(time, list(width, temp), mean)
